home *** CD-ROM | disk | FTP | other *** search
-
-
-
- ____TTTTQQQQLLLL2222((((3333FFFF)))) ____TTTTQQQQLLLL2222((((3333FFFF))))
-
-
-
- NNNNAAAAMMMMEEEE
- TQL2, STQL2 - EISPACK routine. This subroutine finds the eigenvalues
- and eigenvectors of a SYMMETRIC TRIDIAGONAL matrix by the QL method. The
- eigenvectors of a FULL SYMMETRIC matrix can also be found if TRED2 has
- been used to reduce this full matrix to tridiagonal form.
-
-
- SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
- ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
-
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
- rrrreeeeaaaallll dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
-
-
-
- DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
- On Input
-
- NNNNMMMM must be set to the row dimension of two-dimensional array parameters
- as declared in the calling program dimension statement.
-
- NNNN is the order of the matrix.
-
- DDDD contains the diagonal elements of the input matrix.
-
- EEEE contains the subdiagonal elements of the input matrix in its last N-1
- positions. E(1) is arbitrary.
-
- ZZZZ contains the transformation matrix produced in the reduction by TRED2,
- if performed. If the eigenvectors of the tridiagonal matrix are desired,
- Z must contain the identity matrix. On Output
-
- DDDD contains the eigenvalues in ascending order. If an error exit is made,
- the eigenvalues are correct but unordered for indices 1,2,...,IERR-1.
-
- EEEE has been destroyed.
-
- ZZZZ contains orthonormal eigenvectors of the symmetric tridiagonal (or
- full) matrix. If an error exit is made, Z contains the eigenvectors
- associated with the stored eigenvalues.
-
- IIIIEEEERRRRRRRR is set to Zero for normal return, J if the J-th
- eigenvalue has not been
- determined after 30 iterations. Calls PYTHAG(A,B) for sqrt(A**2 +
- B**2). Questions and comments should be directed to B. S. Garbow,
- APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
-
-
-
-
-
-
- PPPPaaaaggggeeee 1111
-
-
-
-